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1. INTRODUCTION 

Presently, it is not well understood how to best model nitrogen oxides (NO x ) emissions from lightning because lightning is highly variable. Peak 
current, channel length, channel altitude, stroke multiplicity, and the number of flashes that occur in a particular region (i.e., flash density) all 
influence the amount of lightning NO x produced. Moreover, these 5 variables are not the same for ground and cloud flashes; e.g., cloud flashes 
normally have lower peak currents, higher altitudes, and higher flash densities than ground flashes [see (Koshak, 2009) for additional details]. 
Because the existing satellite observations of lightning (Fig. 1) from the Lightning Imaging Sensor/Optical Transient Detector (LIS/OTD) do not 
distinguish between ground and cloud flashes, which produce different amounts of NO x , it is very difficult to accurately account for the 
regional/global production of lightning NO x . Hence, the ability to partition the LIS/OTD lightning climatology into separate ground and cloud 
flash distributions would substantially benefit the atmospheric chemistry modeling community. NO x indirectly influences climate because it 
controls the concentration of ozone and hydroxyl radicals in the atmosphere. The importance of lightning-produced NO x is emphasized 
throughout the scientific literature (see for example, Huntrieser et al. 1998). In fact, lightning is the most important NO x source in the upper 
troposphere with a global production rate estimated to vary between 2 and 20 Tg(N)yr' (Lee et al., 1997), with more recent estimates of about 6 
Tg(N)yr' (Martin et al., 2007). In order to make accurate predictions, global chemistry/climate models (as well as regional air quality models) 
must more accurately account for the effects of lightning NO x . In particular, the NASA Goddard Institute for Space Studies (GISS) Model E 
(Schmidt et al., 2005) and the GEOS-CHEM global chemical transport model (Bey et al., 2001) would each benefit from a partitioning of the 
LIS/OTD lightning climatology. 


In this study, we introduce a new technique for retrieving the ground flash fraction in a set of N lightning observed from space and that occur 
within a specific latitude/longitude bin. The method is briefly described and applied to CONUS lightning that have already been partitioned into 
ground and cloud flashes using independent ground-based observations, in order to assess the accuracy of the retrieval method. The retrieval 
errors are encouragingly small. 
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Figure 1. Sample LIS/OTD lightning climatology. 


2. LIGHTNING OPTICAL PROPERTIES 


In the study by Koshak (2007) the frequency distributions of several optical properties of lightning (i.e. flash radiance, flash area, flash duration, # 
of optical groups in a flash, # optical events in a flash) as observed from space were determined for both ground and cloud flashes; it was found 
that the distributions overlapped considerably thereby making it difficult to use the space-based optical measurements to discriminate between 
ground and cloud flashes. However, it was also shown that the means of these distributions are quite different for ground and cloud flashes. 

need to analyze the means has motivated the flash discrimination process discussed in section 3. Fig. 2 shows a sample of the flash area 
distribution for both ground and cloud flashes. We have extended the Koshak (2007) investigation by obtaining the distributions of the maximum 
# events in a group (MNEG) and the maximum group area (MGA). The distributions of MNEG for ground and cloud flashes are shown in Ig. ' . 
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Figure 2. Distributions of flash area. 
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Figure 3. Distributions of MNEG. 
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Figure 4. Sample ground flash showing optical events in space and time. 


Fig. 4 shows an OTD ground flash. The upper right panel is a plan view of 
the optical emissions in the flash, and the lower left panel shows the number 
of optical events within each 2 ms instrument frame. Note that the ground 
flashes often have > 1 frame with > 7 events whereas cloud flashes do not. 


3. THEORY FOR RETRIEVING GROUND FLASH FRACTION (a very brief discussion) 

Consider a set of k - 1, ... , n characteristics given by (x x 9 ..... x . ) for the z th observed flash from space; such a characteristic could 
for example be flash area, or flash radiance, or MNEG. The average of the k th characteristic for N flashes in a latitude/longitude bin is 
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We don’t know N (the # of ground flashes in the bin), but one can obtain solutions in terms of the means (which are estimated using the 
population mean estimates/i , obtained from distributions like those in figures 2 & 3; i.e., we invoke the Central Limit Theorem of statistics): 
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4. APPLICATION, RESULTS, & SUMMARY 

We used National Lightning Detection Network (NLDN) data to partition 5 years of CONUS-only OTD flashes into ground and cloud 
flashes. We then picked 52 locations across CONUS. At each location we picked the radius of a circle that enclosed 1000 OTD flashes (see 
Fig. 5). The circles are colored only to help distinguish the circles for plot clarity. We then applied the retrieval equations in (4). For ri> 1 we 
obtained good results (i.e., small retrieval errors in a ) when characteristics like flash duration and flash # groups were not included. Overall, 
our best results were obtained when we used the simple retrieval equation (n = 1) with the MNEG characteristic . 

in over the 52 regions averaged 0.086 . or just 8.6% of the full range 0-1 of a (see Fig. 5 for specific values of R). Note in figure 5 that 
the value of a shown is the true mean value across the 52 regions; this mean value increases from left (figure 5a) to right (figure 5c) since 
the number of cloud flashes is intentionally depleted from left to right. In summary , the results shown here clearly indicate that reasonable 
retrievals of sround flash fraction can be obtained over CONUS when the population mean estimate of MNEG ( the maximum # of optical 
events in an optical group) is used. Since MNEG is a “peak return stroke detector ” parameter that does not vary sisnificantlv over CONUS , 
then it misht not vary sisnificantlv over the slobe , in which case this retrieval process would also work globally. 


Figure 5a 


Number of Ground Flashes (Red) = 12122 
Number of Cloud Flashes (Blue) = 107166 
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